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^ ■ Abstract 



Using a technique wliich employs a pair of solid scintillator regenerators, the 
E773 collaboration has measured several CP violation parameters in K meson 

decay at Fermilab. We report new results for the phase of r]^ , the Kl — Kg 

mass difference, the Ks lifetime, and the phase difference Arg{r]QQ)—Arg{r]^ ) 

m. K ^ TTTT decay. In addition, we report a measurement of the magnitude 
and phase of r/-| ^ in K tt^tt^j decay. Our preliminary results are com- 
pared with theoretical expectations based on CPT symmetry. 
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I. INTRODUCTION 



Local quantum field theories are automatically invariant under the combined operations 
of charge conjugation, parity, and time reversal: CPT []T|. However, generalizations of 
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quantum mechanics which include gravity might allow pure states to evolve into mixed 
states 0. This sort of unusual (Planck-scale) dynamics could lead to CPT- violating effects 
such as the existence of small differences in particle-antiparticle masses and lifetimes. 
The magnitude of this CPT-noninvariance might be proportional to the ratio of some low 
energy scale to the Planck mass (Ref. [^]). For example, one might have 



for the proton-antiproton mass difference, and 

- ^ = 4.1x10-20 

f^K ^Planck 

for the neutral K system. 

The 90% confidence level limits on the proton-antiproton and electron-positron fractional 
mass differences are approximately 4 x 10"®, twelve orders of magnitude larger than the 
"interesting" region 0. The experimental situation is considerably brighter, however, in 
the sector. Let us define 

AmviKL — > vr+TT") AmviKL vr°7r°) 

= -A 77? V ' ^00 ~ 



Amp(Ks vr+TT ) Amp{Ks tt^tt'^) 

0+_ = Arg{r]+_) , 0oo = Arg{rioo) , A0 = 0oo - 

It can be shown ^ that 

rriK Ktuk / 3 

Here, Am = 0.5286 x W^hs'^c-^ = 3.4793 x IQ-^eV/c^ is the Kl - Ks mass difference. 
The small value for Am provides considerable "leverage" in testing CPT since ^m/niK = 
6.99 X 10^^^. Recently published values yield |m/fo — m-j^\/mK ^ 2.5 x 10"^*^, less than 
two orders of magnitude away from the domain in which Planck-scale physics might play a 
role. 

In Fermilab Experiment 773 we measured 0^ , 0^, and A0 by studying decays of neutral 

K mesons into tt^tt'^ and tt+tt- final states. In addition to our CPT-related tttt results, we also 



report an improved measurement of fJ+^-y, a CP-violation parameter in Kl 7r+7r~7 decays. 
We determined the phases of decay amphtudes by observing the time-dependent interference 
between Ks and Kl decays after a pair of Kl beams passed through regenerators. The 
detector was a reconfigured version of the E731 spectrometer (Ref. 0) which had been used 
to measure Re{e'/e) and other parameters of the neutral K system. 

E773 wrote physics-quality data for about two months, beginning in late July, 1991. 
During this time the experiment recorded ~400 million triggers on nine hundred 8 mm 
cassette data tapes. We present preliminary results from analysis of these data. In Sec. || we 
review relevant neutral K phenomenology. In Sees. |T| and |I^ we describe the experimental 
technique and instrumentation used in E773. In Sees. and ^ we discuss reconstruction of 



candidate K ^ nn and K mi'j events. Simulation of the detector is discussed in Sec. |VII| . 



We describe analysis and fits to our data in Sec. [VIII| . Our conclusions are presented in 
Sec. 1^ These results comprise the thesis work of R.A. Briere and B. Schwingenheuer from 
the University of Chicago, and J.N. Matthews from Rutgers University. 

II. K MESON PHENOMENOLOGY 

A. CP and vrvr decays 

The major source of CP violation in the neutral K system [0,0 is the time-reversal 
asymmetry in the transition rate between K^ and K'^: r{K^ — > K^) > r{K^ — > K^). 
Consequently, the eigenstates of the neutral K system which do not mix in vacuum will 
both contain an excess of K^ relative to K^. We define 

' [2(1 + |e|2)]V2 ' I ^/ [2(1 + |e|2)]i/2 

The phase of e (Ref. [^) is 

(f), ^ tan-^ I — ^1 = (43.33 + 0.14)° , 
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so that Re{e) is positive. The Ks hfetime is ts = 0.8922 x lO'^^s 0. The CP eigenstates 
are 

Neglecting nor mahzat ions, we have \Ki) ~ \Ks) —elKi) and \K2) ~ \Kl) —€\Ks)- A beam 
which is initially pure K2 will evolve into a mixture of K2 and as the Ks component 
of the original K2 state decays. This CP- violating mixing K2 ^ Ki, combined with the 
(CP-allowed) decay Ki rni, provides the dominant contribution to the Kl vttt decay 
amplitude. 

An additional source of CP violation in Kl hit decay is possible if Amp{K2 tttt) is 
nonzero. Since \K2) ~ \K^) — \K'^), this direct CP violation is only possible if Amp{K^ — >■ 
tttt) 7^ A'mp{K'^ tttt). More specifically, we define the isospin and 2 decay amplitudes 
Aq, Aq, A2, A2 so that 

exp^^o Ao = {{7T7t)j=o\T\K^) , exp^^" = ((7r7r),=o|T|:^) 
exp'^'A2 = ((7r7r),=2|T|ir°) , exp'^^ = ((7r7r),=2|T|:^) . 

We find that 

e'^^ exp^(^-^°) ^Al^^ 
V2 {Ao + Ao) 

must be different from zero for there to be direct CP violation in K decay. Here, 60 and 
82 are the (measured) isospin and 2 tttt scattering phase shifts |]12| at the K^ mass. CPT 
invariance requires Aj = Aj*, so the phase of e' is (p^' = 62 — 60 + 90° ~ (47 ±5)° (Ref. [^). It 
may seem surprising that the phase of e' is well known, even though its magnitude remains 
an open question (Ref. [|,|ll|]). By coincidence, the phases of e and e' are nearly equal. If 



we assume CPT invariance, we can write rj^- and r^oo in terms of e and e' as 

r/+„^e + e', r/00 ^ e - 2e'. (2.1) 
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B. CP and iriT'y decays 



CP considerations play a similar role in K ^ 7r+7r^7 decays. The Ks rate is 

dominated by an inner bremsstrahlung amplitude in which one of the final-state pions 
radiates a photon. In Kl decays, the inner bremsstrahlung contribution is CP-suppressed, 
permitting observation of the CP-allowed direct emission process [ri||. As is the case for 



vrvr decays, the CP- violating mixing K2 ^ Ki, combined with the (CP-allowed) decay 



Kl TTTcy, contributes to the Tin'j decay rate in a beam ||15l. We define the ratio of the 
inner bremsstrahlung decay amplitudes: 

^ Amp{KL 7r+7r~7)iB 
^^"^ ~ Amp{Ks 7r+7r"7)iB 

Unless there is an unusually large amount of direct CP violation in Kl 7i~^7i~'j decays, one 
expects /y+.-y ~ ?7+- ~ e. We only observe interference between the inner bremsstrahlung 
amplitudes; the Ks direct emission amplitude is too small to be detected in our data. 

C. CPT 

Several different avenues for CPT violation suggest themselves. One is a CP-violating, T- 
conserving mixing of K^ and K^ to produce the mass eigenstates Kl, Ks- This corresponds 
to a situation where the Ks contains, for example, an excess of K^ while the Kl contains 
an excess of K^. We can describe this by including the CPT violation parameter A in our 
description of Kl, Ks as follows: 

\Ks) ~ (l + e + A)|ir°) + {l-e-A)\W) 
\Kl) ~ (l + e-A)|/CO) - {l-e + A)\W). 

Note the sign switch between e and A in the Kl expression. The component of A which is 
perpendicular to e in the complex plane corresponds to a K^ — K'^ mass difference, while 



the component parallel to e corresponds to a lifetime difference [|TB|. In terms of e, e', and 
A, 
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e + e' - A 



e - 2e' - A . 



(2.2) 



A nonzero value of A will shift r/oo and ri^_ in the same direction in the complex plane. 
From Eq. |2.2| , we find A = e — (2?7_|__ + i]oo)/3. 

Another possible route to CPT violation is through an unusual relationship among the 
various — > vrvr decay amplitudes. Because Aj — Aj* = 2iIm{Ai) and Ai + Ai* = 2Re{Aj), 
the phase of e' would be shifted from its value of (47 ± 5)° by a CPT-disallowed relationship 
such as Aj 7^ Aj* . This would split ?7oo and ?7+_ apart in the complex plane. Other terms 
associated with CPT violation in the decay amplitudes would shift r/oo and ?7_|._ in the same 
direction, much as was the case with A. From Eq. p.2| , we see that e' ^ — i]oo)/3, even 
in the presence of additional (CPT violating) terms which are common to both r^oo and ri^_. 



We extracted phase information about the ri's by measuring the interference between the 
Kl and Ks decay amplitudes in our neutral beam. Since the E773 detector was many Ks 
lifetimes downstream of the production target, we reintroduced a small Ks component into 
our twin beams by passing both beams through regenerators made of plastic scintillator. 

Because the forward elastic scattering amplitudes f{k) and f(k) for and with 
wavenumber k are unequal, a entering a block of material will evolve into a mixture of 
Kl and Ks- The ratio of the Ks and K^ amplitudes a small distance Az inside a regenerator 
will be 



Here, N is the number of scatterers per unit volume. This coherent regeneration amplitude 
p results from constructive interference among all the outgoing spherical waves produced at 
each of the scattering sites in the regenerator. For a finite length regenerator, p includes an 
additional geometric factor associated with the Ks lifetime and the K^ — Ks mass difference. 



D. Ks regeneration 




9 



171 The momentum-dependence of |(/ — has been measured experimentally [T^ and 



is in agreement with the Regge theory prediction that 



' f-T 
k 



oc p"e-*^(2+")/2 . (2.3) 



The connection between a and the phase of / — / comes from analyticity, and is independent 
of Regge theory. 

A kaon which collides with a single nucleus will scatter into a state containing a different 
combination of and K^. Consequently, elastic scattering is regenerative: a pure Kl 
will scatter into a mixed K^, Ks state. Using the optical theorem, one can show that the 
contribution to the Ks amplitude from elastic scattering interferes destructively with the 
scattered kaon's Ks amplitude which arises from coherent regeneration. This cancellation 
is nearly perfect in a regenerator of two interaction lengths. 

We were interested in studying the interference between the (coherently regenerated) 
Ks and Kl decay amplitudes. Since E773's regenerators were shorter than two interaction 
lengths, we made small corrections for contamination in our coherently regenerated K sam- 
ple from "diffraction regeneration" associated with single (and multiple) elastic scattering. 
Regeneration is discussed more fully in the Appendix. 



E. Decay rate downstream of a regenerator 

The Ks and Kl decay amplitudes will interfere downstream of a regenerator. Because 
of the differences in the Kl, Ks masses and decay rates, this interference evolves with proper 
time T. For a state \K) ~ \Kl) + p\Ks), 

where Tl, are the Kl, Ks decay rates. The vrvr decay rate downstream of the regenerator 
is proportional to 

\p\^e-^^^ + |r/pe-^^^ + 2|pr7|e~(^^+^^)^/2 cos(AmcV/^ + 0p - <Pn) ■ 
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The phase of p is approximately -32°, though it varies with energy by several degrees since 
the Kl phase advances with respect to the Ks phase as a kaon travels through a regenerator. 

is 1.12 X 10^°s-\ Tl is 1.93 x lO^s'S and /^mc^/h is 30.67l0-^°s, about 5° per meter 
for 100 GeV kaons. The last term in the decay rate, the interference between Kl and Kg 
decay amplitudes, contains information about 0oo or 4>+-- Defining r_|__ as the K — > 7r"'"7r~ 
decay rate. Too as the K — > 7r°7r° decay rate, and as 0oo — 0+-, one may write 

r+_ - |p|^e-i-^2r^ |r;+_|V-~+ 2|p77+_|e-56^cos(30.6T-77)° 
Too ~ |?7oore-™^+ 2|p?7oo|e-^^"cos(30.6r-77- A0)° 

where r is in units of 10~^° seconds and the arguments of the cosines are in degrees. The 
interference term in Too is most sensitive to changes in when the cosine's argument is 
near -90°, corresponding to r ~ 0. For r ~ 2.5 x 10~^°s the argument of the cosine is near 
zero and the 7r°7r° decay rate is insensitive to small changes in A0. 

III. EXPERIMENTAL TECHNIQUE 

If experimental acceptance and detection efficiency are known with sufficient accuracy, 
A0, and Arg{e) can be determined in a single beam through a measurement of the 

proper time evolution of Too and r+_. However, resolution smearing, diffraction regenera- 
tion, and experimental acceptance influence the reconstruction of charged and neutral final 
states in different ways. As a result, the E773 apparatus was designed conservatively to 
minimize possible systematic errors in the measurement of A0 associated with acceptance 
modeling, energy calibration, and imperfect knowledge of physics parameters such as e'. (We 
subsequently found that our ability to model the behavior of the detector was adequate to 
make single-beam measurements.) 
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A. Double beam technique 



We employed a double-beam technique similar to that used by E731 in its measurement 



a thick regenerator about eleven meters further upstream. This separation was chosen to 
make the vrvr decay rate in the upstream regenerator's beam insensitive to A0 after the beam 
had traveled about a dozen meters downstream of the upstream regenerator, corresponding 
to proper times ~ 2.5 x 10^ The two beams were separated vertically; individual kaons 
passed through only one of the two regenerators. Data were recorded simultaneously for tt^tt^ 
and TT^TT^ decays in both beams. The regenerators switched beams after each accelerator 
spill, approximately once per minute. Various systematic effects associated with detection 
efficiency, kaon flux, and acceptance cancel when rates in the two beams are compared. 



A conservative analysis of the data could be performed by forming a double ratio of the 
decay rates in the two beams for the neutral and charged decay modes. For a given value 
of kaon energy E and decay position z, we deflne 



Here, "up" and "down" refer to the beams containing regenerators in the upstream and 
downstream positions. Binned in E and z, R is insensitive to differences in vr^vr^ and tt+tt" 
acceptance and reconstruction efficiency. For z values downstream of the thin regenerator, 
rQo(z,-E) is nearly independent of A0 for typical kaon energies while V^^^^lzjE) changes 
by about one percent per one degree change in A0. 

A regenerator which is approximately one interaction length in thickness gives the largest 
amount of Kl/Ks interference, since p is proportional to regenerator length, while the 
fraction of beam surviving passage through the regenerator falls exponentially with length. 
Most of the coherent Ks flux in the "upstream" beam, proportional to p^, will have decayed 



(Ref. [T^) of e'. One Ki beam passed through a thin regenerator while the other traversed 



B. Choice of regenerator configuration 
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before traveling as far as the downstream regenerator. A judicious choice for the thickness 
of the downstream regenerator will produce a decay distribution in this beam with similar 
time dependence, except for the proper-time offset in the interference cosine's argument. 
By choosing regenerator lengths which give similar proper time spectra, the importance of 
systematic errors caused by resolution effects is reduced. In particular, the double ratio 
R{z, E) will be close to unity for a wide range of decay positions and energies, and will 
vary from unity by typically 1% per degree of phase difference between 7700 and ?7+-. The 
optimum regenerator thicknesses proved to be 0.4 interaction lengths for the downstream 
and 1.2 interaction lengths for the upstream regenerators. 

A drawback of a double ratio analysis is that it discards events with decay positions 
upstream of the downstream regenerator, sacrificing statistical power in order to reduce 
systematic uncertainties. With adequate modeling of the signal and backgrounds, it is 
possible to extract additional information about the r] phases from events ignored by the 
double ratio technique. We found that we were able to understand (and simulate) the 
behavior of the E773 detector with sufficient accuracy to overcome the systematic effects 
that were of concern at the time the detector was configured. Consequently, we report 
results from the higher precision analysis which includes data from the region between the 
regenerators. 



IV. E773 DETECTOR AND CALIBRATION 



The E773 detector pO[ was designed to determine decay positions and momenta of K 
mesons decaying into 27r, 7r7r7, -kcv, and 37r final states. We used its measurements to 
reject background events in the vrvr sample from semileptonic and 37r decays, and from 
Ks produced through processes other than coherent regeneration. K mesons produced in 
a beryllium target traveled through a series of collimators and sweeping magnets before 
reaching the regenerators. Data from a drift chamber spectrometer and lead glass electro- 
magnetic calorimeter provided kinematic information about charged particles and photons 
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in the final state. Signals from scintillation counters contributed to the experiment's triggers 
and vetoes. A diagram of the detector is shown in Fig. |1]; individual subsystems will be de- 
scribed below. The thin scintillator planes labeled "T,V" inside the vacuum were removed 
midway through the run. The transition radiation detector (TRD) downstream of the last 
drift chamber was not used in the analysis of E773 data. 

A. Neutral beam 

Fermilab delivered a spill of 800 GeV protons to the E773 target about once per minute. 
The protons arrived in 2 ns-long "buckets" spaced by 19ns; each spill lasted twenty seconds 
and contained typically 1.6 x 10^^ protons. A neutral beam was formed by placing sweeping 
magnets, followed by a 5.8m-long copper collimator, downstream of the 36 cm long target. 
The collimator contained two tapered channels which created a clear path for particles 
traveling 4.8 milliradians to the east of, and slightly above or below, the incident proton 
beam's direction. Beryllium blocks placed downstream of the two-hole collimator, together 
with the non-zero targeting angle, reduced the beam's neutron content relative to its kaon 
flux. Some of the beryllium blocks shadowed only the downstream regenerator, reducing 
the hadron flux striking it relative to that in the beam containing the upstream regenerator. 
A 7.6 cm thick lead block (13.6 radiation lengths) converted photons to electron-positron 
pairs which were removed from the beam by subsequent collimators and sweeping magnets. 
With the exception of gaps in the vacuum system at the regenerators, the beams traveled 
in vacuum from 17 m downstream of the target until arriving at the first drift chamber, 
159 m from the target. By the time the beams reached the lead glass array, 181 m from the 
target, each beam was a square, approximately 8 cm on a side. The beams were separated 
vertically with ~15 cm clearance between their inner edges, and passed through holes in the 
lead glass array, stopping in the 3.2 m steel muon filter. Fig. |^ shows the kaon beam profiles 
at the lead glass, determined from Kl tt+tt^ decays. (As with most of the figures, the 
histogram represents data while the superposed points indicate the predictions of the Monte 
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Carlo simulation.) Approximately 10^ K^s (and an equal number of neutrons) struck the 
regenerators each spill. The Ki energy spectrum is shown in Fig. the detector's acceptance 
allowed us to use K mesons with energies above ~ 25 GeV and below ~ 160 GeV. 



B. Regenerators 

In order to allow efficient rejection of events in which a Ks or was produced through 
inelastic processes, both regenerators were built entirely from small blocks of plastic scin- 
tillator. The upstream regenerator was 118 cm long with one quarter of its blocks viewed 
by photomultiplier tubes. The downstream regenerator was 40 cm long with each of its 
blocks viewed by a pair of phototubes. Signals from some of the phototubes participated in 
the experiment trigger; all signals were measured by ADC's and recorded with event data. 
Straight-through muons were used to cross-calibrate the regenerator channels. The moving 
machines which positioned the regenerators (and the beryllium absorber which shadowed 
the downstream regenerator) were controlled by one of the data acquisition computers. The 
regeneration amplitudes for 80 GeV kaons were pup ~ 0.02e~*'^^° and Pdown ~ 0.007e~*'^^° in 
the upstream and downstream regenerators; the magnitude of p decreased with increasing 
kaon momentum as ~ p^^-^. 



C. Tracking system 



The drift chamber spectrometer consisted of two pairs of wire chambers which mea- 



sured particle trajectories upstream and downstream of a large analyzing magnet. The 
chambers' cells were hexagonal, with a maximum drift distance of 6.35 mm. Half-cell offsets 
between paired sense planes provided resolution of "left-right" ambiguities. The chambers 
used an argon-ethane gas mixture which included a small amount of alcohol vapor. Single- 
hit TDC's with 1 ns resolution digitized discriminated chamber signals in response to an 
experiment trigger. Chamber resolution was better than 100 pm per sense plane; frequent 
magnet-off runs of muon data were used to generate alignment constants for the tracking 
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system. Helium bags placed between the chambers served to reduce the effects of multiple 
scattering. Average chamber efficiency was typically 99% per sense plane, but a small num- 
ber of wires showed lower efficiencies. The effects of inefficient wires were modeled in the 
Monte Carlo simulation of the detector. 

The analyzing magnet's vertical field imparted a transverse momentum kick of 
200 MeV/c. Its field integral had been mapped on a two inch grid for E731; the mag- 
net was run at the same field strength in E773. The magnet's aperture was 1.46 m high by 
2.5 m wide by 1 m deep and contained a helium bag during running. 

The tracking system's momentum resolution was typically better than 1%, with a mo- 
mentum dependence given by 

^ = y'(0.45%)2 + (0.012p%)2 

for the momentum p measured in GeV/c. The energy scale was verified by comparing the 
reconstructed masses in X — > 7r"*'7r~ and A — > p7r decays with the known K and A masses. 

D. Lead glass 

The energies and positions of electrons and photons were measured by a lead glass 
calorimeter built from 804 blocks of Schott F-2 glass. Individual blocks were 5.8 cm square 
by 61 cm deep and were stacked parallel to the beam direction in a nearly circular array. 
The radius of the calorimeter was about 0.92m; the beams passed through a pair of holes 
near the center of the array. The array was 18.74 radiation lengths deep. Over the course of 
the run, blocks near the beam holes became less transparent because of radiation damage. 
Periodically, it was necessary to "cure" the blocks with exposure to ultraviolet hght from a 
mercury vapor lamp which resulted in significant (but not complete) recovery. The entire 
calorimeter was located in a temperature controlled, light-tight enclosure. 

An Amperex 2202 photomultiplier tube was mounted on the downstream face of each 
lead glass block to collect Cerenkov hght from electromagnetic showers. Wratten 2A filters 
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between blocks and phototubes transmitted light with wavelengths greater than 430 nm. 
(For shorter wavelengths, the absorption in the lead glass was substantial and varied rapidly 
with wavelength.) Optical fibers, mounted on the upstream faces of all blocks, carried light 
from a single xenon flash lamp which allowed us to monitor detector performance. All blocks 
were constantly illuminated by a small amount of light from a light-emitting diode which 
served to minimize the rate dependence of phototube gains. 

Signals from the lead glass photomultiplier tubes were integrated for 150 ns and digitized 
by LeCroy 1885 ADC's. These dual-range ADC's had conversion gains that were typically 
5 MeV per ADC count below 17 GeV, and 40 MeV per count above 17 GeV. Information 
from channels with >5 counts was recorded with event data. 

Since errors in the measurement of photon energies induced errors in the determination 
of 7r°7r° decay vertex positions, it was necessary to understand in detail the calorimeter's 
behavior. We calibrated the lead glass array with 20 million electrons from the decay K — > 
T:ev (i^^es)- A clean electron sample was obtained after eliminating possible backgrounds 
from A — >• pTT, K — >• tt+tt", and K — >• 7r"'"7r~7r° decays with kinematic cuts. K — >• tt/xi/ {K^^) 
decays were discarded by requiring that no signals be present in the muon veto system. 
Additional cuts required that tracks were well-reconstructed and well-separated at the lead 
glass. Cuts on transverse shower shape and the maximum distance between the cluster 
position and the projected electron track helped eliminate events in which electrons lost 
energy through bremsstrahlung emission. 

The reconstruction of electromagnetic showers depended on three parameters which had 
to be determined for each block: two effective gains and a parameter quantifying a block's 
opacity to Cerenkov light. An iterative procedure was used to determine calibration con- 
stants for each block; the procedure looped repeatedly over electrons in a calibration sample, 
updating constants after each pass, until parameter values stabihzed. The lead glass cali- 
bration parameters changed with time during the run. Central blocks experienced the most 
radiation damage and needed to have their "constants" updated most frequently. Fortu- 
nately, electrons hit blocks near the center of the array most frequently, allowing us 
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to track the most rapidly changing parameters. We were able to generate new calibration 
constants for central blocks from approximately one day's running. 

Fig. 1^ shows the final ratio of the calorimeter energy and track momentum (E/p) for 
electrons. Superposed on the figure is the Monte Carlo simulation's E/p spectrum (see 
Sec. |V11| ); the r.m.s. width of the peak is approximately 3%. The distributions agree except 
for the low side tail in the data which is not yet well-represented by the Monte Carlo. The 
stability of our calibration over time is shown in Fig. ^. The period of time spanned by the 
plot represents about nine weeks, the length of E773's data run. Fig. |^ shows the results 
of fits for the mean and width of the E/p distribution as a function of track momentum. 
Events were restricted to the region near the peak; the apparent shift of the mean away 
from unity is an artifact caused by this restriction. 

E. Trigger and vetoes 

The principal E773 triggers and vetoes were designed to select events with two tracks 
( "charged mode" ) or four photons ( "neutral mode" ) in the final state. Additional triggers, 
useful for calibration and monitoring purposes, selected single muon events, six photon 
events, and "accidentals," events which gave us an indication of the level of additional 
activity in the detector which might be expected to accompany signals produced by the 
decay products of real K mesons. 

Some parts of the veto system participated actively in triggering decisions for both 
charged and neutral mode triggers. For example, signals in the downstream scintillators in 
either regenerator (which were expected to be quiet during coherent regeneration) would 
block both modes' triggers. Other parts of the veto system rejected triggers caused by 
random detector activity as well as triggers from events such as K^-^ decays which were 
accompanied by muons. (We used electrons from the copious decay Kl —>■ neu to calibrate 
the lead glass and to verify the accuracy of the vr+vr" Monte Carlo simulation.) Vetoes for 
the neutral trigger allowed us to reject events in which energy depositions in the lead glass 
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came from hadrons striking the calorimeter. Other neutral mode vetoes permitted us to 
discard triggers in which significant electromagnetic energy went down the beam holes in 
the lead glass array. Analogue signals from most veto counters were measured by ADC's and 
recorded with event data. We used this information during offline analysis to discard events 
in which photons might have missed the lead glass array, based on activity in the photon 
veto system. To reduce backgrounds from inelastic Kg production, neutral and charged 
mode events were discarded when accompanied by activity in either regenerator. 

The charged mode used a two-level trigger. In the first half of the run, the level 1 trigger 
required a signal in any counter in the T,V scintillator hodoscope in coincidence with at 
least two signals from scintillators in the B hodoscope and at least two from counters in the 
C hodoscope (See Fig. 0). The pattern of hits in the B and C counter banks was required 
to be consistent with a pair of tracks from a two-body decay in one of the beams. A signal 
in the /i2 hodoscope (located downstream of a 3.2 m-thick steel wall) would veto a first 
level trigger, allowing us to reject K — > 7r/iz/ decays. The presence of signals in photon 
veto counters immediately before the lead glass would also block the trigger, allowing us to 
discard two-track events which might have come from K tt+tt^tt^ decays. 

The level 2 charged mode trigger in the first half of the run was based on the pattern of 
struck wires in the drift chambers. The trigger required that the four chambers each have 
at least one hit in both left and right sides, as expected in a two-track final state. Midway 
through the run, the T,V hodoscope was removed and a more sophisticated second level 
trigger installed. The new trigger required that the pattern of hits in all four chambers was 
consistent with that expected from a two-body, two-track event, with the possible inclusion 
of extra hits from accidental activity. Both tt+tt^ and 7r+7r^7 final states satisfied the charged 
trigger; no requirement was made on lead glass energy. 

The neutral mode also used a two-level trigger. The first level required that the total 
energy in the lead glass be above ~ 25 GeV. Outputs from groups of nine contiguous lead 
glass blocks were summed by "adders;" individual adder outputs were integrated for 30 ns 
and digitized by ADC's. (This shorter gate helped us eliminate "out of time" clusters 
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registered by the lead glass ADC's.) These sums were combined to form the total energy, 
Et, as a fast analogue sum of signals from the entire lead glass array. The Et > 25 GeV 
requirement was effective since most accidental activity in the lead glass came from muons 
or out-of-time showers, which tended to deposit small amounts of energy in the calorimeter. 

The second level neutral trigger employed a "hardware cluster finder" (HCF) to count 
the number of isolated depositions of energy (clusters) in the lead glass. [^] Events with 
four clusters were accepted as candidate events. 

The level 1 neutral mode vetoes were intended to reject triggers associated with hadronic 
activity in the lead glass or with electromagnetic activity in the vicinity of the calorimeter's 
beam holes. A twenty-one radiation length lead wall behind the lead glass prevented elec- 
tromagnetic showers in the calorimeter from registering in the /il hodoscope. Since hadronic 
showers in the lead glass (which could fake a four-photon trigger) would usually produce 
activity in /xl, a large signal in fil served to veto a level 1 neutral trigger. 

A twenty-eight radiation length lead-lucite shower counter, the BA ("back anti"), was 
installed in the beams behind the calorimeter. Electromagnetic showers tended to deposit 
most of their energy in the upstream part of the BA, while interactions of beam neutrons and 
kaons produced large signals towards the back of the BA. Events with substantial activity in 
the upstream part of the BA, but minimal energy in the last third of the BA, were vetoed. 
The comparison of the signals in the up- and downstream sections allowed us to reject events 
with lost photons without vetoing good events with accidental activity associated with the 
interaction of a beam hadron in the BA. Muons and K^s electrons provided calibration 
information about the BA. 

The CA ( "collar anti" ) covered the inner halves of the blocks around the two beam holes 
and consisted of eight radiation lengths of material followed by scintillation counters. The 
CA allowed us to veto an event in which a photon would have landed near a beam hole; the 
reconstruction of the energy and position of these photons would have been problematic. 

Most of the photon veto counters were constructed of a plane of scintillator followed by a 
lead-lucite shower counter six radiation lengths thick. Each counter's scintillator was viewed 
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by one phototube, while the shower counter was viewed by a pair of phototubes. These 
counters were designed to detect low energy photons which were traveling at angles up to 
50 milliradians with respect to the beams; they were typically 90% efficient at registering 
150 MeV photons. We mapped the position dependence of the response of photon veto 
counters with muons. We were able to calibrate their overall response using K 7r+7r^7r° 



decays in which one photon missed the lead glass. The direction and energy of the missing 
photon could be inferred from the kinematics of the charged pions and reconstructed photon. 
Often this missing photon passed through one of the photon veto counters. By comparing 
the observed signals in the veto counters to the predicted energy of the "undetected" photon, 
the gains and resolutions of the counters were extracted. We included many of the photon 
veto scintillation counters in the level 1 vetoes for charged and neutral triggers. However, 
only the shower counters in the plane closest to the lead glass ("LGA") participated in the 
level 1 vetoes. Information from the shower counters was used during offiine analysis to 
discard events in which photons might have missed the lead glass. 

In addition to the vr+vr" and tt'^tt^ triggers, several special purpose triggers were collected 
at the same time. The "accidental" trigger was used to study the effects of activity in the 
detector not associated with a decaying K meson. Muons traveling through a scintillator 
telescope aimed at the target pile would assert the trigger. The telescope was roughly 50 m 
from the target; the line from the target through the telescope passed approximately ten 
meters west of the spectrometer. As a result, the accidental trigger rate was proportional 
to the instantaneous beam intensity, but was otherwise independent of the presence or 
absence of activity in the spectrometer. Pedestal triggers were used to record the zero- 
signal behavior of all ADC's; calibration flasher triggers were used to monitor the behavior 
of the lead glass calorimeter. Single muon triggers written during special calibration runs 
were used to monitor the performance of many subsystems in the detector. 



21 



F. Data acquisition 



A large number of signals were registered by the detector's data acquisition electronics for 
each event which satisfied the experiment trigger. Discriminated signals from scintillation 
counters were stored in CAMAC latches, while analogue information from the lead glass, 
regenerators, and veto system was processed by FASTBUS ADC's. All data, including 
drift chamber hit times and status information from the trigger electronics, flowed into a 
PANDA [21 1 data acquisition system which could record approximately 13,000 events per 
spill. PANDA wrote data simultaneously to four 8 mm cassette tapes and passed a fraction 
of the events to a micro VaX which monitored the experiment and generated histograms. 
Data describing magnet currents and targeting information were stored at the end of each 
spill along with integrated counting rates from various detectors and components of the 
trigger. 



V. RECONSTRUCTION OF vr+vr- AND vr+vr ^ EVENTS 

A. Tracking 

Tracks in the horizontal (x) and vertical (y) views were found independently. Since the 
analyzing magnet deflected charged particles in the horizontal plane, x-view track segments 
were found separately in the pairs of drift chambers upstream and downstream of the magnet. 
Tracks in the y-view were found in all four chambers simultaneously. Trial track segments in 
the upstream chambers were required to point in the general direction of the decay volume, 
while downstream segments were required to point towards the lead glass. Hits were required 
in at least three of the four sense planes in upstream and downstream x-view track segments. 
A y-view track was required to have hits in at least five of eight sense planes. No tracks 
were allowed to share hits. After track segments were identified, drift distance information 
from the TDCs was used to refine the tracks' positions. 

Track segments in the x-view were projected to the "bend plane" at the center of the 
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analyzing magnet. Upstream and downstream segments were paired if their projections 
were separated by less than 1 cm. The x-view and y-view tracks were correlated using 
information from cluster positions measured in the lead glass. After matching of horizontal 
and vertical views, hit positions were recalculated to correct for small chamber rotations, 
gravitationally-induced sagging of sense wires, and differences in signal propagation time 
along the sense wires. At this point, upstream and downstream y-view track segments were 
refit separately since the analysis magnet tended to impart a small vertical impulse to tracks. 
The difference in upstream and downstream x-view track slopes, combined with knowledge 
of the analyzing magnet's field integral, provided us with sufficient information to calculate 
a particle's momentum. 

Since the paired sense planes in the drift chambers were staggered by a half-cell, the sum 
of the measured drift distances in one view of a chamber should equal the 6.35 mm separation 
between adjacent sense wires. (Because of the 1.1 cm separation between sense planes along 
the beam direction, a small correction was applied for track angle with respect to the z axis.) 
The sums-of-distances allowed us to recognize signals which came from out-of-time tracks. 
Given the temporal structure of the extracted Fermilab proton beam (see Sec. |1V Aj ), a track 



from an earlier (or later) bucket would have its sums-of-distances mismeasured by at least 
1.8 mm. 

We required pairs of tracks to come from a common vertex. The vertex position was de- 
termined using a full fit which included individual tracks' estimated accuracy in contributing 
to the vertex determination. Our requirement on the tracks' distance of closest approach 
depended on track momenta and on the distance from the vertex to the chamber system. 
Typically, tracks with a few millimeter separation satisfied the vertex-quality criterion. Ver- 
tex resolution was ~0.8 mm in the x and y views, and ~ 15 cm in z. The vertex position 
in 7r7r7 decays was set equal to the tracking vertex, without reference to the reconstructed 
photon. 

Additional tracking cuts were made to ensure reconstruction quality. A cut rejected 
tracks with large deviations of hits from ideal positions. A matching ambiguity cut dis- 
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carded events in which the lead glass information could not be used reliably to pair tracks 
with similar x positions with the corresponding y-view tracks. We also required reasonable 
agreement in the projections to the center of the magnet for up- and downstream track 
segments in both x and y views. 

B. Event reconstruction 

Events which came from K tt+tt^ and K — » 7r+7r^7 decays were expected to re- 
construct with invariant mass close to the known kaon mass. Since we were interested in 
detecting kaons which did not scatter in the regenerators, we required that the final state 
TT+TT^ or 7r+7r^7 have minimal transverse momentum (pt) with respect to the incident kaon 
direction. Almost all K^^^ decays were removed at the trigger level so most charged mode 
triggers came from K —>■ ireu decays. There were also significant contributions to the event 
sample from K —>■ tt+tt^tt'' and A —>■ pn^ decays. Events in which a kaon scattered in the 
regenerator contributed to observed backgrounds. 

Backgrounds from Kes decays were suppressed by rejecting events which had a ratio of 
calorimeter energy and track momentum (E/p) close to unity. Requirements on E/p for the 
TT+TT^ and 7r^7r~7 analyses are described in Table |. At the expense of losing a few percent of 
the events with two charged pions in the final state we were able to reject nearly all the K^s 
background. We required tracks to pass through regions of the lead glass in which cluster 
energies could be reconstructed reliably. Events with a track near the outer edges of the 
glass, in a beam hole, or in the collar anti were discarded. 

To reduce background contamination from K^^ decays, we required both tracks to project 
into the fi2 hodoscope. Since "soft" muons were likely to stop in the 3.2 m steel muon filter 
rather than reaching fi2, we rejected events in which either track had momentum below 
7 GeV/c. This cut also reduced the fraction of tt+tt^ decays which were vetoed because a 
pion decayed in flight with the resulting muon striking a /i2 counter. 

Only the most energetic A's from our target survived long enough to decay downstream 
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of the first regenerator. Most of these interacted in a regenerator or were eliminated by 
the requirement that tracks neither hit the collar anti nor went through a beam hole in 
the lead glass. The majority of A —>■ pn decays came from hadronic interactions in the 
regenerators which produced a A without asserting the regenerator vetoes. Because of the 
small Q value, the ratio of the proton and pion momenta in high energy A decay will always 
satisfy Pp/pn > 3.09. Since we had no particle identification for separating hadronic species, 
we assumed that the track with greater momentum was a proton (or antiproton) when 
reconstructing a two-track state as a A (or A). Kinematic cuts (to be described below) were 
used to reject A decays. 

To reduce sensitivity to imperfections in our simulation of accidental activity in the 
detector, both the tt'^tt~ and 7r'*'7r~7 analyses required that reconstructed tracks satisfied 
the charged mode trigger. In addition to this trigger "verification," we discarded events 
where tracks passed near limiting apertures. 

Values for various cuts in the n'^ir' and n^n^'j analyses are listed in Table |. 

C. 7r+7r~ signal and backgrounds 

Reconstructed K vr+vr" mass distributions (after all other cuts, including p"^) are 
shown in Fig. ^. Events from the upstream and downstream regenerator beams are plotted 
separately. In addition, data from the periods before (set 1) and after (set 2) the removal 
of the T,V hodoscope are plotted separately. Results from the Monte Carlo are superposed 
as points on the data histograms. The data's high-side tails contained contributions from 5- 
rays which were not simulated by the Monte Carlo. Our mass resolution was ^ 3.2 MeV/c^ 
in data written before the T,V hodoscope was removed and ~ 0.25 MeV/c^ better after 
its removal; the analysis required 484 MeV/c^ < m^^ < 512 MeV/c^. The tail on the low 
side of the peak came from radiative K ir-K'-f decays. Backgrounds from semileptonic 
decays were nearly negligible. There was no background contribution from tt^tt^ti^ decays 
since the invariant mass of the vt+tt^ system was well below the signal region. Events which 
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had 1110 MeV/c^ < mp^r < 1122 MeV/c^ when reconstructed as A — pvr candidates were 
discarded. 

We calculated the apparent transverse momentum of a kaon with respect to the incident 
beam particle's direction in order to reject backgrounds from diffractive and inelastic K 
production as well as semileptonic decays. We extrapolated the kaon trajectory from the 
decay vertex upstream to the regenerators. The kaon's transverse momentum, p^, is the 
component of the observed vector momentum which is perpendicular to the line connecting 
the production target with the kaon's position at a regenerator. Because of the spacing 
between the beams' inner edges we could determine unambiguously which regenerator the 
kaon traversed. 

The distributions are shown in Fig. The vrvr analysis required < 250 (MeK/c)^; 
tails near the coherent regeneration peak in the distribution were dominated by radiative 
K TTTT'y decays and diffractively regenerated Kg- At large p^, inelastically produced Ks 
contributed significantly to the tails. Some smearing of the data's coherent peak was due to 
the presence of 5-rays. 

The p'^ distributions of events passing all other cuts were used to determine background 
levels. The small backgrounds consisted of a mix of diffractive K tt+tt^ decays, radiative 
7r+7r~7 decays, and semileptonic final states; they are described in Table |T|. A total of 1.8 
million tt~^tt~ decays were found in the entire run's data after all cuts; the event sample is 
described in Table We required that kaons had energies in the range 30 to 160 GeV 
and that decay vertices were downstream of the vacuum window following the appropriate 
regenerator. To reduce backgrounds associated with interactions of beam hadrons in material 
near the downstream regenerator, we discarded events from the upstream regenerator beam 
with decay vertices near the gap in the vacuum system. Figs. ^ and |10] show the decay vertex 
z position and K energy spectrum for 7r"'"7r~ decays passing all other cuts. Data (histogram) 
and Monte Carlo (points) are shown for each beam. The T,V hodoscope, located near 
z = 141 m, was included in the charged mode trigger for set 1 data. 
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D. Tr~^ IT 7 signal and backgrounds 



The 7r"''7r~7 analysis required two well-reconstructed tracks and at least one in-time lead 
glass cluster with energy >1.5 GeV which was not associated with either track. In the center 
of mass of the kaon, this photon was required to have an energy greater than 20 MeV. To 
suppress background from tt+tt^tt^ decays, we required 

, 2 



T \2 



= pI < -0.01125 



4[(p^_)2 + m^_] 

Here, m+_ is the invariant mass of the two charged tracks and is the square of 

their transverse momentum with respect to the parent kaon. (We assumed that the kaon 
traveled in a straight line from the production target to the decay point and did not scatter 
in the regenerator.) Since po is proportional to the vr" momentum component in the kaon 
center-of-mass which is parallel to the beam direction, Pq> for events which might be Sir 
decays. [^,^ Shown in Fig. ^ are energy spectra for the photon in the K center of mass 
for 7r+7r^7 decays satisfying all other cuts. 

We removed background from K vr+vr" events accompanied by accidental activity in 
the calorimeter by discarding events with m+_ > 484MeV/c^. Remaining backgrounds from 
A —>■ piT decays (in combination with an accidental lead glass cluster) were reduced with 
kinematic cuts. Events which reconstructed to have 1100 MeV/c^ < < 1130 MeV/c^, 
total energy above 100 GeV, and momentum ratio Pp/pn > 3.0 were discarded. 

To select events in which only coherently regenerated Ks and decay amplitudes 
participated, we rejected events with reconstructed mass outside the range 484 MeV/c^ < 
m7r+7r-^ < 512 MeV/c^ or with large transverse momentum such that p"^ > 150 (MeV/c)^. 
Shown in Figs. |1^ and ^ are mass and for ir^n'j decays satisfying all other cuts. Results 
of the Monte Carlo simulation (shown as points) did not include background modeling. 

The data separated naturally into two sets corresponding to the periods of running time 
before and after the removal of the T,V hodoscope. After cuts and background subtractions, 
a total of (10,769 ± 17) 7r+7r~7 decays were found in the two data sets; the distribution of 
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events between the sets is shown in Table |V. 

Events for each set were binned according to kaon momentum p and decay position z 
in a grid with cell size 10 GeV/c x 2 m. We required kaons to have energies between 25 



and 155 GeV. Figs. and show the decay vertex z position and K energy spectrum 
for 7r+7r~7 decays passing all other cuts. Decay vertices were required to be downstream of 
the vacuum window following the appropriate regenerator. The mass and distributions 
of events passing all other cuts were used to determine the level of background in each 
set. Backgrounds were typically 2%, and are listed for each subset in Table 0. The small 
backgrounds consisted of a mix of noncoherent K 7r+7r^7 decays and final states made 
from two-track events combined with accidental activity in the calorimeter. We modeled the 
shape of the background using events which passed the mi'j mass cut, but had between 
300 and 2000 (MeV/c)^. These events were binned in {p,z) grids identical to those of the 
signal region, scaled appropriately, then subtracted from the signal-region {p, z) grids. 



VI. RECONSTRUCTION OF -k^t:^ EVENTS 
A. Photon reconstruction 

The reconstruction of an electromagnetic shower ( "cluster" ) in the lead glass array pro- 
ceeded in stages. Lead glass channels whose ADC signals were local maxima were taken 
as the central blocks in candidate clusters. Signals from central blocks were required to be 
above 200 ADC counts; signals from the eight blocks surrounding each central block were 
included in the energy estimate formed for clusters. 

A first estimate of a cluster's energy was made using the sum of the ADC counts, divided 
by the effective gains of each channel, in the cluster's nine blocks. The effective gain included 
photomultiplier gain, ADC digitization conversion gains, photocathode quantum efficiency, 
and light collection efficiency. A rough correction for absorption of Cerenkov light in each 
block was applied to the initial cluster energy estimate. With this correction, the estimated 
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cluster energy was usually within 30% of the true shower energy. 

Several corrections were made to this initial estimate. A small amount of energy was 
added to account for blocks which did not meet the readout threshold of five ADC counts. 
A second correction compensated for radial leakage outside the nine blocks in a cluster, into 
the beam holes, or past the outer edges of the array. When showers overlapped, the energy in 
shared blocks was divided between showers based on typical shower profiles. During the early 
part of data taking, the lead glass enclosure's temperature control system malfunctioned. 
Since we recorded the enclosure's temperature with event data, we were able to correlate 
gain changes with temperature, and to map the enclosure's temperature as a function of 
time. Using this information, we applied a small time-dependent gain correction to lead 
glass signals. The ADC pedestals showed a small rate-dependent shift; we corrected them 
based on the instantaneous event rate. 

The final determination of the energy in an electromagnetic shower included correction 
for effects associated with absorption of Cerenkov hght in the lead glass. The depth at 
which a shower deposited most of its light increased with shower energy. As a result, 
the response of the calorimeter tended to "brighten" with increasing shower energy, since 
less hght was absorbed by the glass between "shower max" and the photomultipher tube. 
The absorption coefficients varied from block to block, and changed with time as radiation 
damage increased. The coefficient of the cluster's central block was used to calculate the 
correction factor. An additional correction compensated for the fact that photons began 
showering deeper in the lead glass than did electrons. Shower energies associated with 
photons which converted to e+e" pairs in the B and C hodoscopes were adjusted by 3.5%. 
In a final correction we shifted electron-induced shower energies 0.3% to compensate for the 
effects of bremsstrahlung photons produced in the trigger hodoscopes near the calorimeter. 

Shower x, y positions were found using the estimated energies deposited in the three 
columns and three rows of blocks in a cluster. We studied the calorimeter's position reso- 
lution by comparing the impact point of electrons, determined by the drift chambers, with 
shower centers, determined by the lead glass. The position algorithm's results were typically 
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accurate to 3 mm, though the resolution varied from 1.5 mm near the corners of blocks to 
5 mm at the centers of blocks. 



B. Event reconstruction 



The invariant mass of two photons originating from a single decay vertex and traveling 
at small angles relative to the beam direction is 



where z is the distance from the decay vertex to the calorimeter, Ei and E2 are the photons' 
energies, and is their separation at the calorimeter. We were interested in reconstructing 
final states in which two (or three) 7r°'s produced four (or six) detected photons. There were 
three possible ways to "pair" the photons in a candidate 7r°7r° event, and fifteen ways in a 
candidate 37r° event. For each of the possible pairings, we calculated a separate decay z for 
each pion by assuming that the invariant mass of the corresponding photon pair equaled the 
7r° rest mass. We then determined errors for each z from the uncertainties in cluster energies 
and positions, and calculated a ior the hypothesis that the photon pairs came from the 
same vertex. Using the pairing with the smallest x^, we assigned a decay z for the event 
based on a weighted average of the pions' individual z positions. Typical vertex resolution 
was about 1 meter; the invariant mass of the four (or six) photons was 



Several cuts were made to reduce background contamination in the data. We rejected 
events with pairing > 4. In addition, we discarded four-photon events in which the 
second-best pairing had a reasonable ^^'^ yielded an invariant mass near the K mass. 
Background from 37r° decays with undetected photons was reduced by discarding events with 
significant activity in the photon veto system. In addition, 37r° decays with photon energies 
below the HCF thresholds were removed by cutting events with > 600 MeV of in-time energy 
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in any lead glass block which was not associated with one of the four reconstructed clusters. 
Backgrounds associated with interactions of beam hadrons with material in the detector 
were reduced by cuts on the number of drift chamber hits and on trigger hodoscope activity. 
As was the case with charged mode events, we rejected tt^tt^ candidates with significant 
activity in either regenerator. 

Other cuts were made to simplify Monte Carlo simulation of the detector. Photons were 
required to have energies greater than 2.2 GeV to reduce sensitivity to inaccuracies in the 
simulation of HCF thresholds. Events were discarded if any cluster reconstructed within a 
half block of the outer edge of the lead glass array. We rejected events with appreciable 
activity in the collar anti. To avoid problems from events with overlapping photons, we 
required the energy distribution among the blocks in a cluster to be consistent with the 
shape expected for a single electromagnetic shower. 

The calorimeter was calibrated with electrons, not photons. The coupling between the 
energy and z scales in the invariant mass relationship [Eq. ( |6.1| )] was used to check the photon 
energy scale. The reconstructed positions of the downstream faces of the regenerators were 
studied as functions of energy; the energy scale was adjusted so that the data and Monte 
Carlo distributions agreed. 



Cuts made in the vr vr and A0 analyses are summarized in Table VI. 



C. vr vr signal and backgrounds 

Our preliminary result for Argiri^^o) — Argij]^ ) is based on data written after the T,V 

hodoscope was removed. (See Sec. [IV E| .) The reconstructed K n^-n^ mass distributions 
(after all other cuts) are shown in Fig. |16|. Results from the Monte Carlo simulation of 
signal and backgrounds are superposed on the data. In the plots on the left, the simulation 
includes both signal and backgrounds. In the plots on the right, only the simulated back- 
grounds are superposed on the data. The signal region is between 474 MeV and 522 Mev. 
The background is mostly from K Stt*^ with two undetected photons. There is a small 
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background contribution from neutron interactions in detector material which produced two 
neutral pions. 

The X, y vertex position was not measured in neutral decays so the transverse momentum 
of the parent kaon, pt, could not be determined directly. Instead, the center-of-energy of the 
four photons was used to select decays of kaons which had not scattered in the regenerators. 
The center of energy was defined as 

where fi was the point at which the i^^ photon struck the calorimeter. Since fc.E. corre- 
sponded to the point at which the kaon would have reached the lead glass if it had not 
decayed, unscattered kaons had tc.e. inside the beam spots. 

We defined an event's "ring number" to be the area circumscribed by a square, centered 
on the closer beam, whose perimeter contained fc.E. (see Fig. |1^). The ring number distri- 
butions (after all other cuts) are shown in Fig. |I8[ Unlike the charged mode spectra, the 
neutral mode ring plots contain no sharp peaks because of the finite beam size. Events were 
required to have a ring number less than 120 cm^. Near the coherent regeneration region in 
the ring plot is a shoulder from diffractively regenerated Ks', inelastically regenerated Ks 
populate the high ring number region. As before, the predictions of the Monte Carlo are 
superposed on the data. 

Normalizations for our background subtractions were established using events which lay 
outside the signal region in the mass vs. ring number plane. The background shape in the 
signal region was known from Monte Carlo studies of ^ Svr^ decays which faked 7r°7r° 
events and from measurements of the p"^ distribution for tt+tt^ events. (Our p"^ resolution 
in the charged decay mode let us identify scattered kaons whose centers-of-energy remained 
inside the beam profile.) A total of 2.6 x 10^ tt'^tt" decays were found in the set 2 data after 
all cuts. The event sample is described in Table |V11| . The kaon energy was required to be 
in the range 40 < Ex < 150 and the decay vertex z position was required to be in the 
range 120 < z < 152 for the upstream regenerator and 130 < 2; < 152 for the downstream 
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regenerator. The vertex position distributions and kaon energy spectra for tt^tt^ decays 
passing all other cuts are shown in Figs. |19| and |20. 



VII. MONTE CARLO SIMULATION 

We measured the phases of ?7+_, ?7+_^, and 7700 by studying the proper time evolution 
of the TTTT and inc'-f decay rates downstream of the regenerators. The detector's acceptance 
varied with kaon energy and decay position, so we needed to unfold apparatus effects from 
observed decay spectra in order to determine values for the magnitudes and phases of the 
?7's. We used a detailed Monte Carlo simulation of the experiment to calculate the p and z 
dependence of the detection efficiency. We also used the Monte Carlo to study backgrounds 
and to determine the effects of various sources of systematic uncertainty. The accuracy of 
the simulation could be verified with the large samples of Kes, tt+tt^tt^, and Svr^ decays 
which satisfied the experiment triggers. 

The Monte Carlo program generated raw data which was processed by the same recon- 
struction programs used to analyze real data. The effects of spurious detector activity could 
be studied by combining data from accidental triggers (see Sec. |1V h^ ) with simulated event 
data. 



A. Beamline simulation 

The neutral beam simulation modeled the creation of and in our production 
target as well as the propagation of kaons through absorbers and the (imperfectly aligned) 
collimation system. Even the small amount of residual Ks amplitude present in high energy 
kaons striking our regenerators was modeled by the beam simulation. Propagation of kaons 
in the regenerators was treated in detail, allowing for coherent regeneration, diffraction 
regeneration through single (and multiple) elastic scattering, inelastic processes, and decays 
inside the regenerator. Regeneration in air gaps, vacuum windows, and the T,V hodoscope 
was included in the Monte Carlo. The simulation of the interactions of daughter particles 
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from K decay included multiple Coulomb scattering, the production of bremsstrahlung 
photons by electrons and positrons, and the conversion of photons to e^e~ pairs. Simulated 
charged pions were allowed to decay in flight while traveling through the spectrometer. 
We based our , K'^ energy spectra on Malensek's parametrization of the and K~ 



production spectrum for 400 GeV protons incident on a beryllium target ||2^ . He described 
the differential production cross section in terms of xp = PK/Ebeam and pt', in general, 
cross sections expressed in terms of these variables scale well with beam energy (Ref. 0). 



Because of measurement errors in the result of Ref. and uncertainties associated with 
the extrapolation to higher incident beam energy, we found it necessary to adjust the shape 
of the K production spectrum to fit our data. The correction function varied from unity by 
~ ±20% over the range of usable kaon energies. Though it was very convenient to have the 
Monte Carlo's energy spectrum well-matched to the data's energy spectrum, our analysis 
was not particularly sensitive to inaccuracies in spectrum modeling. 

The Monte Carlo propagated kaons through the target and the beryllium and lead ab- 
sorbers (see Sec. [IV A| ), modeling the effects of coherent regeneration and (non-regenerative) 
single elastic scattering in addition to attenuation of the beam. The simulation modeled the 
X, y intensity profile of the incident proton beam and allowed kaons to scatter in collimator 



walls. The accuracy of the absorber/collimator description can be seen in Fig. ^ for Kf,3 
decays. Tails in the data's distribution of angle between the target-to-vertex direction and 
the z axis are well-matched by the Monte Carlo. These tails are associated with elastic 
scattering of kaons in beamline elements. 

The X, y locations and sizes of limiting apertures in the detector were measured using 
electrons from Kes decays, while the z positions of the apertures were determined using 
survey information. Shown in Fig. ^ are comparisons of the data and Monte Carlo electron 
illuminations at the east edge of the HDRA aperture in 1 mm bins. Illumination of charged 



pions at the HDRA and lead glass are shown in Figs. ^ and ^ for events which passed 
TT+TT^ analysis cuts. 
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B. Drift chamber simulation 



The Monte Carlo simulation of the drift chambers did not attempt to model detailed 
processes like electron drift and diffusion. Instead, the drift time distribution observed in 
the data was integrated to produce a mapping of drift distance -^-^ drift time. Chamber 
resolution was incorporated by jittering the true positions of particles at chamber sense 
planes with the appropriate Gaussian distributions. Smeared positions were then converted 
to drift times. When two tracks passed through the same drift cell, only the first signal 
to arrive at the chamber amplifier/discriminator was included with simulated data. (Our 
chamber TDC's only registered the first hit received by each channel.) 

Chamber inefficiencies were modeled in some detail, since they determined the charac- 
teristics of the chamber component of the charged mode trigger. Individual wire efficiencies 
were ~99%; the characteristics of especially inefficient channels observed in the data were 
included in the simulation. 

The chamber simulation did not model S-recy production. Consequently, the high-side 
tail in the track distribution is not simulated. We expected that the effects of 5-rays 
were not correlated significantly with the decay vertex z position, and would only cause an 
overall loss of data relative to Monte Carlo without introducing biases in our phase-of-77 
measurements. 

C. Lead glass simulation 

The mean amount of Cerenkov light which reached an individual block's photomultiplier 
tube was a function of the incident particle's energy, the block's attenuation length for 
Cerenkov light, and the depth at which the shower began. Since showers were almost fully 
contained in the calorimeter, signals from showers which began deep in a block suffered less 
attenuation than those which began early. We parametrized the distribution about the mean 
for the light reaching the phototube as a Gaussian with high-side and low-side exponential 
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tails. Photon showers were simulated as pairs of electron showers which begin where the 
photon converted to an e^e~ pair. The use of an effective block length allowed us to take into 
account the fact photons began to shower at different depths in the calorimeter. Extensive 
studies using EGS4 have shown this to be a sensible description of the production of 
Cerenkov light in lead glass. ||2^ 

Studies with the xenon calibration flasher were used to determine the mean number of 
photoelectrons per Cerenkov photon expected for each block. The simulation chose the 
amount of Cerenkov light which reached each phototube, determined the expected number 
of photoelectrons, then varied this assuming a Gaussian distribution whose width was the 
square root of the expected mean. 

The event simulation modeled cluster shapes using a library of electron clusters obtained 
from special runs. In these calibration runs, a beam of photons was converted to e+e^ pairs; 
the e"*" and e~ were steered by calibration magnets (see Fig. P into the lead glass. (By 
varying the currents in the magnets we were able to illuminate the entire calorimeter with a 
wide range of electron energies.) Clusters were selected from the library based on incident 
electron energies and positions in the central block. The signal in each block was scaled by 
a constant to give the simulated response. For photons, two appropriate electron clusters 
were overlayed. 

We adjusted the tail structure in the glass response to correct for effects not modeled in 
the EGS4 simulation. For example, EGS4 assumed that the attenuation of Cerenkov light 
was uniform along the length of a block. In reality, absorption increased with depth: there 
was more radiation damage deep in blocks due to hadronic interactions. In addition, the 
curing lamp (which was placed near the upstream face of the calorimeter) was less effective 
at correcting damage in the downstream end of blocks. 
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D. Veto and hodoscope simulation 

Efficiencies of scintillation counters were measured using muon tracks. Small gaps be- 
tween adjacent counters in hodoscope planes were mapped using the drift chambers and 
included in the Monte Carlo. Since we had determined the response of photon veto counters 
with data we were able to simulate their response to photons in the Monte Carlo. 



E. Comparisons of data and Monte Carlo 

The large samples of — > ttcz/ and K Stt'' decays provided us with useful tools for 
verifying the accuracy of the Monte Carlo simulation. We collected roughly five times as 
many 37r° decays as 27r° decays, and more than ten times as many -f^es's as tt+tt^ decays. 
It was particularly important that we modeled correctly acceptance and resolution effects 



that influenced the reconstructed vertex z distributions. Shown in Figs. ^ and are 
comparisons of these distributions for data and Monte Carlo. Our success in simulating the 
acceptance of these "high statistics" modes gives us confidence in the validity of conclusions 
we draw from characteristics of the vrvr and tttt'') samples. 



VIII. ANALYSIS AND FITS 

A. (f)^ , Am, Ts analysis 

We used the Monte Carlo detector simulation to predict the tt+tt" acceptance in 10 GeV/c 
X 2 m cells in momentum p and decay position z. The acceptance in each cell was defined 
as the ratio of the number of detected events with reconstructed p, z in the cell divided by 
the number of generated events with true p, z in that cell. Since the simulation accurately 
modeled the data, effects due to resolution smearing were properly included in our fits. 
Average acceptance was approximately 20%. 
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We measured Arg{'ri^_) as well as the regeneration parameters a and — ^^jto g v 
with a fit to background-subtracted data which included free parameters describing the 
beam intensities and energy spectra. In this fit we constrained the Kl — Ks mass difference 
to the value measured by E731: Am = 0.5286 x 10^^hs^^c~^ (Ref. 0), while using the 
world-average value for the Ks lifetime: ts = 0.8922 x 10~^°s (Ref. [0]). (These values 
correspond to 0^ = 43.33°.) We found that = (43.35 ± 0.70)° where the error is 
statistical (systematic uncertainties will be discussed below). The fit's 598 for 577 

degrees of freedom. 

In a separate fit we determined Am and ts (as well as the beam and regeneration 
parameters) by constraining the phase of 1]+^ to equal the phase of e: 

0+_ = 0^ = tan~^ (2Amc^Ts/fij . 

We found that Am = (0.5286 ± 0.0029) x lO^o^s-^c-^ and Tg = (0.8929 ± 0.0014) x IQ-^'^s, 

corresponding to (f)^ = 43.35°. Ag ain, the errors are statistical. When 0_| was uncon- 

strained, we found that 0+_ - 0, = (-0.84 ± 1.42)°, Am = (0.5268 ± 0.0041) x lO^^fis-^c-'^, 
and Ts = (0.8942 ± 0.0026) x IQ-^^s. 

Sources of systematic errors included uncertainties in our determination of experimental 
acceptance, errors associated with background subtractions, uncertainties in our descriptions 
of the two regenerators, and possible deviations of the momentum dependence of (^f — /k 
from a power law. We studied acceptance issues by moving the edges of limiting apertures 
and introducing spurious 2;— dependent slopes (some of which were functions of momentum) 
to force disagreement between data and Monte Carlo. Background systematics were investi- 
gated by varying the cut, thus changing the amount of K^.^ contamination remaining with 
the coherent signal. We changed the positions and relative lengths of the regenerators in the 
apparatus simulation to determine our sensitivity to uncertainties in regenerator geometry. 
Issues concerning an unexpected momentum dependence in the regeneration amplitude are 
described in Ref. [|rB| . The sensitivity of 04._ to these effects is summarized in Table [Vllij ; 



combining the errors in quadrature yielded a net systematic uncertainty of ±0.79°. The fit 
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for did not allow Am or ts to vary. Increasing the values of Am and ts by one (Par- 
ticle Data Group) standard deviation (Ref. ^) would change 0+_ by +0.38° and —0.62° 
respectively. 

In summary, our preliminary results are 

= [43.35 ± 0.70 (stat.) ± 0.79 {syst.)Y 
(Am fixed at 0.5286 x lO^^hs'^c'^, ts fixed at 0.8922 x 10"^°s.) 
5(0+_) = +0.38° for 5(Am) = 0.0024 
5(0+_) = -0.62° for S{ts) = 0.0020 . 

If instead we constrain 0+_ = tan^^i^lAmc^Ts/h), we find 

Am = [0.5286 ± 0.0029 (stat.) ± 0.0022 (syst.)] x 10^°^s"^c~^ 
Ts = [0.8929 ± 0.0014 (stat.) ± 0.0014 (syst.)] x lO'^^s . 

Fitting simultaneously for Am, and ts yields 

0+_ - 0e = [-0.84 ± 1.42 (stat.) ± 1.22 {syst.)]° 

Am = [0.5268 ± 0.0041 (stat.) ± 0.0029 (syst.)] x lO^^hs^^c-^ 

Ts = [0.8942 + 0.0026 (stat.) +0.0018 (syst.)] x 10"^°s . 

Systematic uncertainties will decrease with additional analysis work. 

These results can be compared with the Particle Data Group's world average values 

1^ and previous measurements from E731 [^ and NA31 [§. The PDG values are (p^ = 

(46.5 + 1.2)°, Am = (0.5351 + 0.0024) x lO^^hs-^c'^, ts = (0.8922 + 0.0020) x lO^^^s. 

Since the value of (p^ is correlated with the values of Am and ts, the PDG used the 

world-average values of Am and Ts when combining 0^ results from different experiments. 

If we were to recalculate our 0^ result using the PDG values for Am and ts we would 

obtain 0+_ = (44.38 + 0.70 + 0.79)°. The E731 results were 0+_ = (42.2 + 1.5)°, Am = 
(0.5286 + 0.0028) x lO^^hs-^c-'^, ts = (0.8929 + 0.0016) x lO^^^s. The resuh for was 

39 



obtained from fits which assumed ts = 0.8922 x 10 but allowed Am to float. The results 
for Am and fixed = (p^- If we were to recalculate our (f)^_ result using the Am value 



obtained [|1|] in E731's 0+_ fit, we would find 0+_ = (42.14±0.70±0.79)°. The NA31 result 
was determined assuming the PDG values ts = 0.8922 x 10^^'^ and Am = 0.5351 x 10^°. 
They found 0+_ = (46.9 ± 2.2)°. This can be compared with our result recalculated 
using the PDG values for Am and ts mentioned above. Our results are consistent with the 
absence of CPT violation in Kl tt^tt^ decays. 

B. A0 analysis 

Our preliminary result for 0oo ~ 4'+- is based on data written after the T,V hodoscope 
was removed, approximately 70% of the full data set. As was the case in the 0+_ analysis, 
we binned data and Monte Carlo samples for both tt^tt^ and tt+tt^ decays in 10 GeV/c x 
2 m (p, z) cells. The copious Stc^ sample served to verify the accuracy of the neutral mode 
Monte Carlo. Acceptance for neutral mode decays varied with position and K energy, and 
was typically ~ 5%. 

We fit our background-subtracted data for values of A0, the regeneration parameters 



a and 



{f ~ jQ Q ^'5 beam intensity and energy spectrum parameters. 



70 GeV 

2 



For the cuts described in Table 0, the fit x was 682 for 619 degrees of freedom. We 
found that A(f) returned by the fit varied somewhat with the value of the minimum cluster 
energy cut, min- For our preliminary result we report the average of two fits which used 
Ed min = 2.2 GeV and E^i min = 4.0 GeV, and assign a systematic uncertainty of 0.5° to 
this effect. We find A0 = (0.67 ± 0.85)°, where the error is statistical. 

Other systematic uncertainties arose from possible inaccuracies in our determination of 
acceptance, our background subtractions, and our parametrization of the lead glass calorime- 
ter. Effects involving regeneration influenced both 0oo and and canceled when extract- 
ing their difference. As was the case with the charged mode analysis, we introduced spurious 
z— dependent slopes to force disagreement between data and Monte Carlo in order to gauge 
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the consequences of acceptance errors. We varied the amounts of background from Svr^ 
decays, inelastically and diffractively regenerated Ks-, and from beam interactions in mate- 
rial in the detector. We studied several possible sources of calorimeter misparametrization, 
including an energy scale mismatch between the lead glass and the tracking system, mis- 
representation of the lead glass' energy resolution, and uncorrected pedestal shifts. The 
dependence of on these effects is summarized in Table |X|; the quadrature sum of the 
various sources of systematic uncertainty is ±1.1°. Consequently, our preliminary result for 
the phase difference between 7700 and 77+ _ is 

A0= [0.67 ± 0.85 {stat.) ± 1.1 {syst.)]° . 

We expect the statistical error to decrease when the full data set is used in the analysis; the 
systematic error will decrease as the analysis continues. 

We can compare our result with the Particle Data Group's world average value 
and previous measurements from E731 and NA31 [Q. The 1992 PDG value is A(/) = 
(—0.1 ± 2.0)°. The E731 result of A0 = (—1.6 ± 1.2)° was determined with a fit which al- 
lowed the values for Am and to vary. If E731 had calculated their result using our Am and 
Ts values they would have obtained A0 = —1.75°. The NA31 result was A0 = (0.2 ± 2.9)°, 
and was determined assuming the PDG values ts = 0.8922 x 10^^*^ and Am = 0.5351 x 10^°. 
Their value is insensitive to changes in Am. Our result is consistent with previous measure- 
ments and with the absence of CPT violation in — > tttt decays. 

C. r/^ ^ analysis 

The Monte Carlo simulation of the apparatus was used to determine the acceptance as 
a function of p and z for TTTc-f decays in each beam. We generated simulated data sets which 
were approximately twenty times larger than the corresponding real data sets and determined 
average acceptance in each cell of the (p, z) grids described above. Overall acceptance of our 
apparatus for the 7r+7r~7 decay mode was typically 11%. The momentum spectrum and flux 
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of kaons striking the regenerators were determined from the more copious K vr+vr" decay 
mode. Characteristics of the regeneration amphtude p were also determined with tt+tt" 
data. The 7r+7r^7 analysis (and Monte Carlo) assumed that Am = 0.5286 x 10^°/is~"'^c~^ 
(Ref. 0]) and that the Kg n+n-'j branching ratio (Ref. [g) was 4.87 x lO'l 

We extracted phase information about rj^^^ from the interference between the inner 
bremsstrahlung Kl^s t^t^I decay amplitudes. The direct emission amplitude (which 
tended to produce a more energetic photon in the kaon rest frame) only participated in the 
Kl decay rate, and did not interfere with the Ks decay amplitude at an observable level. 
Based on a previous measurement (Ref. IQ)) the Monte Carlo simulation assumed that the 
direct emission contribution to the decay rate with E* > 20 MeV was 68.5%. Here, E* 
is the photon energy in the K rest frame. 

The four background-subtracted (p, z) data grids were fit simultaneously for the mag- 
nitude and phase of f]+-j, the Ks 7r+7r^7 branching ratio, and the ratio of the di- 
rect emission and inner bremsstrahlung contributions to the decay amplitude using a 
log-likelihood procedure. Our preliminary result is l^^+.-yl = (2.414 ± 0.065) x 10^^ and 
= Arg{ri^_^) = (45.47 ±3.61)°, where the errors are statistical. After fitting, a was 
calculated using all bins with at least five events. We found = 366.9 for 347 degrees of 
freedom. Results from the individual data sets and the combined sample are summarized in 
Table 0. 

Possible sources of systematic error included uncertainties associated with the regener- 
ation amplitude, background subtractions, and our parametrization of the kaon flux and 
momentum spectrum. In addition, our descriptions of Kl niT'-f direct emission and the 
Ks — * TTTC'-f branching ratio contributed to the systematic error. By varying the analysis pa- 
rameters by one standard deviation, we studied the extent to which systematic uncertainties 
influenced our conclusions. The results are presented in Table Adding systematic errors 
in quadrature yields the preliminary results 

\ri+_^\ = [2.414 ± 0.065 (stat.) ± 0.062 (syst.)] x 10"^ 
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= [45.47 ± 3.61 (stat.) ± 2.40 isyst.)]° . 

These errors are significantly smaller than those of the best previous measurement (Ref. |15| ) 
which found \ri+^^\ = (2.15 ± 0.26 ± 0.17) x lO^^ and = (72 ± 23 ± 17)°. Our value 

for ri^_^, and the best previous measurement, are compared in Fig.p7|. This measurement 

of f]^ ^ is also not significantly different from the analogous tt'^tt" amplitude ratio = 

(2.268 ± 0.023) x 10"^ (Ref. 0) and the 0+_ value presented in Sec. |VIIIA| . This result is 
consistent with the absence of CP violation beyond that present in Kl tttt decays. 

IX. CONCLUSIONS 

We report preliminary results from analysis of data from the 1991 fixed target run at 

Fermilab. Our results for Arg{ri^ ), Am, ts, \ri^_^\, and Arglt]^ ^) are based on the entire 

data set. Our result for Acf) = Arg^rjQo) — Arg{ri+_) is based on analysis of about 70% 
of the available data. Systematic uncertainties will decrease as we refine the analyses; the 
statistical precision of our A(f) result will improve when the full data set is used. 

Our two-pion results are consistent with CPT conservation in K ^ tttt decays. Our 
?7+_-y result is consistent with the absence of an unusual source of CP violation in K ^ nnj 
decays. 

APPENDIX: REGENERATION 

1. Ks coherent regeneration 

Neglecting decays, the wave function for a relativistic meson with wavenumber k will 
obey the Klein-Gordon equation (V^ + k'^)il)K{^ = 0. If the passes through a scattering 
medium, its total wavefunction will become the sum of the incident and scattered waves. 
Defining f(k) to be the average forward scattering amplitude at each scattering site inside 
the medium, we find that ipx must satisfy 
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{V^ + +47iNf{k))ijK = 

with approximate solution 

^^(x,t)^e^[(^+^)'=^-^*/^] (Al) 

for a plane wave traveling in the z direction. Here, is the number of scatterers per unit 
volume. Similarly, the wavefunction for a K'^ with forward scattering amplitude f{k) will 
be 

^(x,t) ^ e L J . (A2) 

Note that these expressions include the effects of elastic and inelastic scattering since a 
positive imaginary part to f{k) (or f{k)) causes lip] to decrease. This is sensible: the optical 
theorem guarantees Imf > whenever the total scattering cross section is nonzero. 

We can use Eqs. (|A1| , |A^) to describe the position dependence of a mixed state traveling 
in the z direction through a regenerator. The 2;-dependence of the Ks, Kl amplitudes for 
the state I A') = As\Ks) + A l\Kl) will obey 



OA 



s 



dz 

BAl 
dz 



ik + iNn 
ik + iNTT 



' f + T 
k 

' f + T 
k 



As + zNttI^^^^Al 
Al + zN7r(^^)As 



neglecting weak interactions and decays. A state which is pure Kl at 2; = will evolve into 
a mixture of Kl and Ks as long as (/ — /) 7^ 0. In particular, the ratio of the Ks and Kl 
amplitudes a small distance Az inside a regenerator will be 

This coherent regeneration amplitude p results from constructive interference among all the 
outgoing spherical waves produced at each of the scattering sites in the regenerator. Since 
the wavelength of a 100 GeV/c K meson is ~ 2 x 10^^^ meters, coherent scattering in a 
macroscopic regenerator only takes place for extremely small scattering angles. 
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The total cross section at high energy is dominated by inelastic processes so the scattering 



amplitudes / and / are nearly imaginary The momentum-dependence of \{f — f)/k 
has been measured |]18| and is in agreement with the Regge theory prediction ||T9[ that 



k 



(A4) 



The connection between a and the phase of / — / comes from analyticity, and does not 
depend on the validity of Regge theory. 



2. Diffraction regeneration 

Scattering does not mix ^ K^. However, since the scattering probability for a 
differs from that for a K^, a K meson which contains a mix of K*^ and K'^ amplitudes will 
scatter into a state containing a different combination of and K^. In terms of the elastic 
scattering amplitudes fk{9,(f)) and ff.{9,(j)), the state 

iJi = A\K°) +A\K^) 

will scatter elastically into the state 



^|Jf = AMe,<f>)\K') +Af,{e,<p)\K^) 



For convenience, we define f22{0, (p) = fk{9, 0) + /fc(6', (p) and /2i(6', (p) = fk{9, (p) - f^{9, 0). 
If we rewrite the incident and final states tpi, ipf in terms of Ks and Kl, we have 

= As\Ks) +Al\Kl) (A5) 

{Asf22 + ALf2l)\Ks) + {ALf22 + Asf2l)\KL)\^. (A6) 

In the final state Ks amplitude, the term Asf22 corresponds to elastic Ks Ks scattering 
while the term A^f 21 comes from (regenerative) K^ Ks scattering. Similarly, in the final 
Kl amplitude, the term A^f 22 arises from elastic K^ — > Ki scattering, while As f 21 results 
from Ks —>■ Kl scattering. 
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A pure Kl will scatter into a mixed Kl, Ks state; we can define the diffraction regener- 
ation amplitude po for elastic scattering into ^,0 as the ratio /2i//22- With this definition 
we can rewrite Eq. ( |A6| ) as 

[As + PdAl)\Ks) + {AL + PDAs)\KL)\f^2^— . (A7) 

Because fk and f^ are nearly imaginary (Ref. [^), /22 ~ ^1/22! so that 

to A\ •^^i . /21 

J22 I /22 I 

With the help of the optical theorem, this can be rewritten as 

p.(».*)«^(^) (A8) 

when 9,(f) are close to the forward direction so that f22{0,(f)) ~ /22(0,0). Here, ar is the 
i^L-nucleus total scattering cross section. 

There are interference effects associated with the presence of both coherent and diffractive 
processes. Imagine that a enters a regenerator of length L, scatters elastically after 
traveling a distance Az, then leaves the regenerator without rescattering. (We will ignore 
effects associated with finite regenerator length such as the decay of the Ks amplitude and 
the phase drift between the kaon's Ks and Kl components due to the non-zero Kl — Ks 
mass difference.) From Eq. (|A3|) , the kaon's wavefunction immediately before scattering will 
be 

l/jbefore = I ^ ll^s) + (A9) 

neglecting normalization. The Ks component is the result of coherent regeneration. Imme- 
diately after scattering, the kaon's wavefunction will become 

j f tN7TAzf2iiO,0) , \ , tN7TAzf2iiO,0) V^.l.e^'^ 

Rafter = < I ^ + pD I \Ks} + I 1 + ^ pD I \Kl} > f22— ■ (AlO) 

The first term in the Ks amplitude corresponds to Ks — > Ks scattering while pl, corresponds 
to elastic Kl — >■ Ks scattering. For the E773 regenerator geometry and K energy range, 

46 



both contributions to the Ks amphtude of Eq. ( |A10|) are smalL The modification to the 



Kl amphtude is second-order in the regeneration parameters /21, po- As the kaon travels 
through the remaining distance {L — /\z) in the regenerator, its Kg and Kl amphtudes 
wiU continue to evolve because of Ks,l Kl s and Kl^s ~^ Kl^s forward scattering. In 
particular, its Ks amplitude will gain an additional contribution from coherent regeneration. 
Neglecting normalization, and terms which are second order in the regeneration parameters, 
we find that the kaon's wavefunction at the downstream end of the regenerator is 

^ = ]( 1 + Pd\\Ks) + \Kl)> f 22^ ■ (All) 



With the definition of the coherent regeneration amplitude p in Eq. (|A3|) , we can rewrite 
Eq. ( [ATT| ) as 

^/ = [{p + Pd)\Ks) + |irL)]/22^. (A12) 

Because of the relative minus sign between p and pr){0,4'), there will be destructive 
interference between coherent and diffractive regeneration. If we approximate /2i(6',0) ~ 
/2i(0, 0) and make use of Eqs. ( [A3D and (|A8|) , we find 

For a two interaction length regenerator, NLaT = 2, and the destructive interference will 
be nearly perfect. 

The interference between the coherent and diffractive regeneration amplitudes is influ- 
enced by a number of factors in addition to regenerator length. If fk{6,(f)) and fi^{6,(f)) 
are not purely imaginary, Eq. (|A8|) will be inaccurate and the cancellation will be im- 
perfect. To the extent that dependence on momentum transfer spoils the approximation 
/22(^5 0) ~ /22(0,0), Eq. ( |A13| ) will be unreliable. The nonzero K^-Ks mass difference 
makes the Ks amplitude's phase lag behind the Kl phase as the kaon travels through the 
regenerator: the Ks amplitude from elastic scatters at the upstream and downstream ends of 
the regenerator will contribute with slightly different phases to the decay amplitudes. This 
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phase drift is negligible at Fermilab energies, where it is ~ 5° per meter of flight path for a 
100 GeV kaon. In addition, the finite Ks lifetime plays a role. The contribution to the Ks 
amplitude from coherent forward scattering will decay in flight as the kaon moves through 
the regenerator. The fraction which survives is a function of the regenerator length and 
the kaon's lifetime, and is independent of the location of the elastic scatter. The diffractive 
contribution to the Ks amplitude, created at the scattering site, will also decay in flight. As 
a result, the surviving diffractive contribution depends on the distance from the scattering 
site to the downstream end of the regenerator. Consequently, the sum of the coherent and 
diffractive Ks amplitudes will vary slightly with the location of the scattering site. 
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TABLES 




TABLE I. Requirements on 


events used in tt+tt" and 7r'*"7r~7 


analyses. Set 1 (2) data were 


recorded before (after) removal of the T,V hodoscope. 






+ — 




Invariant mass, MeV/c 


484 < mTTTT < 512 


484 < mTTTTT < 512 


p|, (MeV/c)2 


< 250 


< 150 


track momentum, (jreV/c 


> 7 


> 7 


E/p 


< 0.8 


< 0.85 


Kaon energy (GeV) 


30 < Ek < 160 


25 < < 155 


Decay vertex, upstr. reg. (m) 


118.5 < z < 127 (set 1) 


117.33 <z< 139.33 (set 1) 




.UK. Iz9 < z < 14U 






lio.o S ^ S J--^' (.set zj 


117.33 <z< 159.33 (set 2) 




.OR. 129 < z < 154 




Decay vertex, dnstr. reg. (m) 


i^n-r" -/ <' 1/in ('00+ 1 ^ 
iou ^ z 2^ i4u (^ser i j 


128.63 < z < 140.63 (set 1) 




130 < z < 154 (set 2) 


128.63 <z< 158.63 (set 2) 


A suppression 


mp^ < 1110 MeV/c^ 


mp^ < 1100 MeV/c^ 




.OR. mp^ > 1122 MeV/c^ 


.OR. mpTT > 1130 MeV/c^ 






.OR. £;a < 100 GeV 






.OR. ^>p/p7r < 3 


TABLE IL Estimate of backgrounds in tt+tt 


sample. 




set 1 


set 2 


Upstream Regenerator 


0.29% 


0.25% 


Downstream Regenerator 


0.90% 


0.75% 
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TABLE III. 


Number of coherent tt^tt events. 






set 1 set 2 


Whole Run 


Upstream Regenerator 


707k 750k 


1,457k 


Downstream Regenerator 


160k 207k 


367k 


Both Regenerators 


867k 957k 


1,824k 


TABLE IV. Number of coherent tt^tt 7 events after background subtraction. 




set 1 set 2 


Whole Run 


Upstream Regenerator 


4,559 ± 11 3,779 ± 9 


8,338 lb 14 


Downstream Regenerator 


1,209 ± 6 1,222 ± 6 


2,431 ± 9 


Both Regenerators 


5,768 ± 13 5,001 ± 11 


10,769 ± 17 


TABLE V. Estimate of backgrounds in tt+tt 7 sample. 




set 1 


set 2 


Upstream Regenerator 


75.5 ± 11.0 events 


65.0 ± 9.1 events 




1.6% 


1.7% 


Downstream Regenerator 


28.0 ± 6.0 events 


31.9 ± 6.3 events 




2.3% 


2.5% 
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TABLE VI. Requirements on events used in the analysis. 













Invariant mass (MeV/c^) 


474 < TO^^ < 522 






484 < m^TTTT < 512 


p|, (MeV/c)2 


— 






< 250 


ring number (cm^) 


< 120 






— 


track momentum (GeV/c) 


— 






> 7 


cluster energy (GeV) 


^cluster ^2.2 






— 


E/p 


— 






< 0.8 


Kaon energy (GeV) 


A<d<EK< 150 






SO < Ek < 160 


Decay vertex (m) 


upstr. reg.: 120 < 


z < 152 




upstr. reg.: 118 < z < 127 










.OR. 129 < z < 152 




dnstr. reg.: 130 < 


z < 152 




dnstr. reg.: 129 < z < 152 


A suppression 








rup^ < 1110 MeV/c^ 










.OR. rup^ > 1122 MeV/c^ 


TABLE VII. Number of coherent 7r°7r° events. 


(Note that only data set 2 is used in the 


preliminary result for Ac/)). 












set 1 




set 2 


Whole Run 


Upstream Regenerator 


85k 




180k 


265k 


Downstream Regenerator 


35k 




75k 


110k 


Both RogencM'ators 


120k 




255k 


375k 
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TABLE VIII. Estimates of systematic error in tt+tt analysis. 



Source of uncertainty change in (p^ 

Acceptance calculation 0.57° 
Analyticity constraint on phase of p 0.5° 
Background subtractions 0.2° 
Regenerator descriptions 0.05° 
Total 0.79° 

Increase Am by 0.0024 (1 PDG a) 0+_ increases +0.38° 

Increase ts by 0.0020 (1 PDG a) 0+_ decreases -0.62° 



TABLE IX. Estimates of systematic error in analysis. 



Source of uncertainty change in Acp 

Acceptance parametrization for 7r''"7r~ and 7r°7r° 0.3° 

Background subtractions 0.6° 

Calorimeter energy scale 0.6° 

Calorimeter nonlinearity (including minimum cluster energy cut) 0.5° 

Calorimeter resolution 0.5° 

Total 1.1° 



TABLE X. ?7-| -y magnitude and phase from different data sets. 



Regenerator 


set 1 




set 2 




Upstream 




= (2.358 ± 0.086) x 10"^ 


l^+-7l 


= (2.406 ± 0.084) X 10-^ 




Arg{r]^ 


__^) = (51.1 ±5.4)° 




__^) = (40.9 ± 6.5)° 


Downstream 




= (2.376 ± 0.098) x 10"^ 




= (2.443 ± 0.083) x 10"^ 






__^) = (43.1 ±7.8)° 




.-y) = (54.0 ±6.5)° 


All data 




= (2.414 ± 0.065 ± 0.062) x IQ--'^ 




__. ) = (45.5 ±3.6 ±2.4)° 
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TABLE XI. Estimates of systematic error in 7r"'"7r 7 analysis. 



Source of uncertainty 


fractional change 


change in r/_| ^ 


change in in (pj^ 


Kg — 7r7r7 ranching ratio 


2.3% 


0.030 X 10-3 


1.76° 


DE/IB ratio 


6.0% 


0.041 X 10-3 


0.19° 


regeneration power law a 


2.5 % 


0.001 X 10-3 


0.92° 


Normalization 


1.6% 


0.021 X 10-3 


1.22° 


Background normalization 


16% 


0.008 X 10-3 


0.15° 


Background shape 




0.030 X 10-3 


0.17° 


Data/MC p spectrum mismatch 


2 X 10"^ 


0.002 X 10-3 


0.54° 


Total 




0.062 X 10-3 


2.40° 
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FIGURES 

FIG. 1. Elevation view of the E773 detector. Kaons in the beams travel to the right in the 

figure. Details of the individual hardware systems appear in the text. The thin hodoscope labeled 
"T,V" was removed partway through the run. 

FIG. 2. Neutral beam profiles at the lead glass array, as determined with Ki —> tt+tt" decays. 
Data are shown as the histogram with error bars; results from the Monte Carlo simulation are 
shown as points. 

FIG. 3. Predicted A^l energy spectrum striking the regenerators. The smooth curve is based 
on the parametrization of and K~ production by Malensek, modified to agree with our data. 

FIG. 4. Ratio of the calorimeter energy and track momentum (E/p) for K^^ calibration elec- 
trons. The r.m.s. width of the peak is approximately 3%. Data (histogram) and Monte Carlo 
(points) are shown. The low side tail in the data (which is not well-represented by the Monte 
Carlo simulation) is under investigation. 

FIG. 5. Mean E/p during the run. Shown in the plot is the average value of E/p for K^s 
electrons as a function of run number. The period of time spanned by the run was about nine 
weeks. 

FIG. 6. E/p mean and r.m.s. width vs. electron track momentum. The upper figure shows 

the results of fits for the mean E/p, restricted to the region near the peak, for Kes electrons. The 
apparent shift away from unity is an artifact of the restriction in the fits to the peak region. The 
lower figure shows the r.m.s. width, again restricted to the region near the peak. Data (histogram) 
and Monte Carlo (points) are shown. 
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FIG. 7. 7r"'"7r~ invariant mass distributions for events satisfying all other analysis cuts. Events 
from decays in the upstream and downstream regenerator beams are plotted separately. Midway 
through the run the T,V hodoscope was removed; data are plotted separately from the periods 
before (set 1) and after (set 2) its removal. Data (histogram) and Monte Carlo (points) are shown. 
The data's high-side tails contain contributions from 5-rays which are not simulated by the Monte 
Carlo. 

FIG. 8. p\ distributions for K vr^vr^ decays satisfying all other analysis cuts. Events from 
decays in the upstream and downstream regenerator beams are plotted separately. Midway through 
the run the T,V hodoscope was removed; data are plotted separately from the periods before (set 
1) and after (set 2) its removal. Data (histogram) and Monte Carlo (points) are shown. Some 
smearing of the data's coherent peak is due to the presence of (5-rays which are not simulated by 
the Monte Carlo. 

FIG. 9. Decay vertex z distributions for K — tt+tt" decays satisfying all other Arg{r]-\ ) 

analysis cuts. To reduce backgrounds associated with interactions of beam hadrons in material 
near the downstream regenerator, events with vertex z between 127 m and 129 m were discarded. 
Data (histogram) and Monte Carlo (points) are shown for each beam. The T,V hodoscope, located 
near z = 141 m, was included in the charged mode trigger for set 1 data. 

FIG. 10. Reconstructed kaon energy distributions for K — ^ 7r"^7r~ decays satisfying all other 
analysis cuts. Data (histogram) and Monte Carlo (points) are shown for each beam in the two 
data sets. 

FIG. 11. Photon energy spectrum in the K rest frame for 7r"'"7r~7 decays satisfying all other 
cuts. Events are required to have E* > 0.02 GeV. Data (histogram) and Monte Carlo (points) are 
shown for the entire data sample. The vertical dashed line indicates the location of the 0.02 GeV 
cut. 
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FIG. 12. 7r''"7r~7 invariant mass distributions for events satisfying all other cuts. Data (shown 
as histograms) are plotted separately for each beam. Results of the Monte Carlo simulation (shown 
as points) do not include background modeling. 

FIG. 13. pj^ distributions for K — > 7r"'"7r~7 decays satisfying all other cuts. Data (histogram) 
and Monte Carlo (points) are shown, plotted separately for each beam. The Monte Carlo did not 
simulate backgrounds. 

FIG. 14. Decay vertex z distributions for K — ^ Tr'^Tr~j decays satisfying all other cuts. Data 
(histogram) and Monte Carlo (points) are shown, plotted separately for each beam in the two data 
sets. 

FIG. 15. K energy distributions for K —>■ 'k~^'k~j decays satisfying all other cuts. Data (his- 
togram) and Monte Carlo (points) are shown. 

FIG. 16. 7r°7r'^ invariant mass distributions for events satisfying all other cuts. Data (histogram) 
and Monte Carlo (points) arc shown. In the plots on the left side of the figure, the simulation 
includes both signal and backgrounds. In the plots on the right, only the simulated backgrounds 
are included in the Monte Carlo sample. 

FIG. 17. Illustration of "ring number" fov a, K ^ tt'^tt^ decay. In this example, an event's 
center of energy in the calorimeter is at (x, y) = (14, 4) cm, closer to the upper beam which is 
centered at {x,y) = (0, 11.6) cm. The square whose perimeter contains the center of energy has 
side length 2 x 14 cm and area 784 cm^. As a result, the event's ring number is 784 cm^. 

FIG. 18. TT^TT^ ring distributions for events satisfying all other cuts. Data (histogram) and 

Monte Carlo (points) are shown. In the plots on the left side of the figure, the simulation includes 
both signal and backgrounds. In the plots on the right, only the simulated backgrounds are included 
in the Monte Carlo sample. 
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FIG. 19. Decay vertex z distributions for K — tttt events satisfying all other cuts. Data 
(histogram) and Monte Carlo (points) are shown. 

FIG. 20. Reconstructed kaon energy distributions for K 7r°7r° decays satisfying all other 
cuts. Data (histogram) and Monte Carlo (points) are shown. 

FIG. 21. Angle between the target-to-vertex direction and the z axis in K^s decays, projected 
onto the horizontal plane. Wings of the distribution correspond to K mesons which scattered in 
beamline elements before decaying. Data (histogram) and Monte Carlo (points) are shown. 

FIG. 22. Comparison of the data and Monte Carlo Kg^ electron illuminations at the east edge 
of the HDRA aperture in 1 mm bins. Data (histogram) and Monte Carlo (points) are shown. 

FIG. 23. TT^ track illuminations at the HDRA. Data (histogram) and Monte Carlo (points) are 
shown. 

FIG. 24. TT^ track illuminations at the lead glass array. Data (histogram) and Monte Carlo 
(points) are shown. 

FIG. 25. Decay vertex z distributions for K Stt^ decays. Data (histogram) and Monte 

Carlo (points) are shown. The large Stt*^ sample provides a valuable check of our vr'^Tr'' acceptance 
calculation. 

FIG. 26. Decay vertex z distributions for K ireu decays. Data (histogram) and Monte Carlo 
(points) are shown. The large i^es sample allows us to check our tt+tt" acceptance calculation. 

FIG. 27. Magnitude and phase of ry-i -y from this experiment and from E731. The smaller error 

bars indicate statistical uncertainties. The larger errors correspond to a sum, in quadrature, of 
the quoted statistical and systematic errors. For comparison, the vertical dashed lines indicate the 
±1(7 bounds on our (f)j^- measurement while the horizontal dashed lines indicate the iblcr bounds 
on the Particle Data Group's world-average value for |?7^ |. 
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